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Ultra high energy cosmic rays: clustering, GUT scale and neutrino masses 
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The clustering of ultra high energy (above 5 • 10 19 eV) cosmic rays (UHECR) suggests that they might be 
emitted by compact sources. We present a statistical analysis on the source density based on the multiplicities. 
The propagation of UHECR protons is studied in detail. The UHECR spectrum is consistent with the decay of 
GUT scale particles and/or with the Z-burst. The predicted GUT mass is mx = 10 6 GeV, where b = 14.6U 7 - 
Our neutrino mass prediction depends on the origin of the power part of the spectrum: m„ = 2.75 /Iq g® eV for 
halo and O. 26 IU 4 eV for extragalactic (EG) origin. 


1. Introduction 

The interaction of protons with the microwave 
background predicts a drop in the cosmic ray flux 
above the GZK jij cutoff «5T0 19 eV. The data 
shows no such drop. Section 2 studies the propa¬ 
gation and determines the probability P(r, E, E c ) 
that protons created at distance r with energy E 
reach earth above a threshold E c . Using this P 
one can give the observed spectrum by one nu¬ 
merical integration for any injection spectrum. 

It is an interesting phenomenon that the 
UHECR events are clustered. Recently, a sta¬ 
tistical analysis J|] based on the multiplicities of 
the clustered events estimated the source density. 
In Section 3 we extend the above analysis. 

In Section 4 we study the scenario that the 
UHECRs are coming from decaying superheavy 
particles (SP) and we determine their masses mx 
by an analysis of the observed UHECR spectrum. 

Ultrahigh energy neutrinos (UHEz/) scatter on 
relic neutrinos (R^) producing Z bosons, which 
can decay hadronically (Z-burst) |J. In Section 5 
we compare the predicted proton spectrum with 
observations and determine the mass of the heav¬ 
iest Rzz via a maximum likelihood analysis. 

The details of the presented results and a more 
complete reference list can be found in 10- 

2. Propagation of UHECR protons 

Using pion production as the dominant effect 
of energy loss for protons at energies >10 19 eV 


ref. p] calculated P(r,E,E c ) for three thresh¬ 
old energies. We extended the results of [jfj. 
In our Monte-Carlo approach protons are prop¬ 
agated in small steps (10 kpc), and after each 
step the probabilities of pion production and the 
energy losses due to pair production, pion pro¬ 
duction and the adiabatic expansion are calcu¬ 
lated. We used the following type of parametriza- 
tion P(r, E, E c ) = exp [—a • (r/1 Mpc) 6 ]. Fig. [I| 
shows a(E/E c ) and b(E/E c ) for a range of three 
orders of magnitude and for five different E c . 

3. Density of sources 

The arrival directions of the UHECRs mea¬ 
sured by experiments show some peculiar clus¬ 
tering: some events are grouped within ~ 3°, 
the typical angular resolution of an experiment. 
Above 4 • 10 19 eV 92 cosmic ray events were de¬ 
tected, including 7 doublets and 2 triplets. Above 
10 2 ° eV one doublet out of 14 events were found 
||. The chance probability of such a clustering 
from uniform distribution is rather small [^]|]]. 

The clustered features of the events initiated an 
interesting statistical analysis assuming compact 
UHECR sources The authors found a large 
number, ~ 400 for the number of sources within 
the GZK sphere. We generalize their analysis. 
The most probable value for the source density is 
really large; however, the statistical significance 
of this result is rather weak. 

Fig. | shows the CL regions for one of our 
models (injected energy distribution c{E) oc P -3 ; 
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Figure 1. The parametrization of P(r, E, E c ). 



Figure 2. The la (68%) and 2cr (95%) con¬ 
fidence level (CL) regions for j* and the source 
density (14 UHECR with one doublet). 


Schechter’s luminosity distribution: h{j)dj oc 
07J* ) ” 1 ' 25 ex P ( (j/j*))- The regions are 
deformed, thin ellipse-like objects. For this model 
our final answer for the density is • 

10 -3 Mpc~ 3 , where the first errors indicate the 
68%, the second ones in the parenthesis the 95% 
CLs, respectively. The choice of [p] h(j)oc <5(j)— 
and, e.g. E~ 2 energy distribution gives much 
smaller value: 2.7772 6 5 3(2 1 7o)10 -3 Mpc -3 , which 
is in a quite good agreement with their result. 

4. Decay of GUT scale particles 

An interesting idea discussed by refs. 00 is 
that SPs could be the source of UHECRs. The 
hadronic decay of SPs yields protons. They are 
characterized by the fragmentation function (FF) 
D(x,Q 2 ) which gives the number of produced 
protons with momentum fraction x at energy 
scale Q. For the proton’s FF at present accel¬ 
erator energies we use ref. |Ii|. We evolve the 
FFs in ordinary and in supersymmetric QCD to 
the energies of the SPs. This result can be com¬ 
bined with the prediction of the MLLA technique, 
which gives the initial spectrum of UHECRs at 
the energy rrix (cf. Fig.. |^). Similar results are 
obtained by [Q. 

Depending on the location of the source -halo 
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Figure 3. The quark FFs at Q=10 16 GeV 
for proton/pion in SM (solid/dotted line) and in 
MSSM (dashed/dashed-dotted line). We change 
from logarithmic scale to linear at x = 0.01. 



Log (E [eV]) 

Figure 4. UHECR data with their error bars 
and the best fit from a decaying SP. There are no 
events above 3 x 10 20 eV (shown by an arrow). 
Zero event does not mean zero flux, but an upper 
bound for the flux, thus, the experimental flux is 
in the ’’hatched” region with 68% CL. 


or extragalactic (EG)- and the model SM or 
MSSM- we study four different scenarios. In the 
EG case protons loose some fraction of their en¬ 
ergies, described by P(r, E, E c ). We compare the 
predicted and the observed spectrums by a max¬ 
imum likelihood analysis. This analysis gives the 
mass of the SP and the error on it. 

Fig. ^ shows the measured UHECR spectrum 
and the best fit, which is obtained in the EG- 
MSSM scenario. The goodnesses of the fits for the 
halo models are far worse. The SM and MSSM 
cases do not differ significantly. The most impor¬ 
tant message is that the masses of the best fits 
(EG cases) are compatible within the error bars 
with the MSSM gauge coupling unification GUT 
scale: nix = 10 b GeV, where b = 14. 

5. Z-burst scenario 

Already in the early 80’s there were discussions 
that the UHE^ spectrum could have absorption 
dips at energies around -E£® s =M§/(2 to^)= 4.2 • 
10 21 (1 eV /m v .) eV due to resonant annihilation 
with Rj/s of mass m v , predicted by the hot Big 
Bang, into Z bosons of mass Mz |l^JI^| . Recently 
it was realized that the same annihilation mech¬ 
anism gives a possible solution to the GZK prob¬ 
lem It was argued that the UHECRs above 


the GZK cutoff are from these Z-bursts. 

We compare this scenario with observations. 

The density distribution of Rz'S as hot DM fol¬ 
lows the total mass distribution; however, it is 
less clustered. Thus we, as opposed to practically 
all previous authors ||[n]|uj, do not follow the 
unnatural assumption of having a relative over¬ 
density of 10 2 -j-10 4 in our neighborhood (for an 
approach with lepton asymmetry see Jl8|] ). 

We give the energy distribution of the produced 
protons in our lab system, which is obtained by 
Lorentz transforming the CM collider results. 

The next ingredient is the propagation of the 
protons, which can be described by P(r, E p , E). 

Finally, we compare the predicted and observed 
spectrum and extract the mass of the Ri/ and the 
necessary UHEz/ flux by a maximum likelihood 
analysis. In the Z-burst scenario small R^ mass 
needs large ££ es to produce a Z. Large E r * s results 
in a large Lorentz boost, thus large proton en¬ 
ergy. In this way the detected energy determines 
the mass of the Rzx The analysis is completely 
analogous to that of the previous section. 

Our best fits to the observed data gives for the 
neutrino mass 2.75+g gi^'ggj eV for the “halo”- 
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Figure 5. Differential v+D fluxes (averaged over 
the families) required by the Z-burst hypothesis. 
The horizontal errors indicate the uncertainty of 
the mass determination and the vertical errors in¬ 
clude also the uncertainty of the Hubble expan¬ 
sion rate. Also shown are upper limits from Fly’s 
Eye and the GLUE, as well as projected sensitivi¬ 
ties of AMANDA, Auger and OWL. 


and 0.263*10 ^ 4(0 22) ^ or “EG”-case, respec¬ 
tively (in the halo/EG cases the power part of the 
UHECR spectrum is generated in the halo/EG). 
This gives an absolute lower bound on the mass 
of the heaviest v of 0.06 eV at the 95% CL. 

The most attractive pattern for v masses is hi¬ 
erarchical. Using the mass difference of the at¬ 
mospheric v oscillation for the heaviest mass |20| , 
one obtains values between 0.03 and 0.09 eV. It is 
an intriguing feature of our result that the smaller 
one of the predicted masses is compatible on the 
ss 1.3 ct level with this scenario. 

The necessary UHEz/ flux can be obtained from 
our fits. We summarized them in Fig. ||, together 
with some upper limits. We also compared our 
7- flux with the EGRET bound ( 2 ^] . Our 7 flux 
is somewhat smaller than that of EGRET. 
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